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Abstract 

We have modified the sexual Penna model by introducing the fluctuating environ- 
ment and fluctuations representing physiological functions of individuals. Additionally, 
we have introduced the mother care corresponding to the protection against the delete- 
rious influence of the environment, the learning capacity of individuals corresponding to 
their immunity and adaptation to the environment fluctuations and the other risk factors 
appearing at puberty. Each of the above mentioned elements influences mainly the sur- 
vival of newborns and young individuals while genetic defects accumulated in the genomes 
increase the noise of individuals and were responsible for higher mortality of older indi- 
viduals. All these modifications enable precise fitting the age structure of the simulated 
populations to the age distributions of human populations. 

1 Introduction 

In 1952 Medawar formulated the theory of ageing called the "mutations accumulation theory" 
PQ|. He argued that: if there is a random death killing the individuals in a population indepen- 
dently of their age, the sizes of groups of individuals at a given age decrease exponentially with 
age; if all individuals keep the same reproduction potential per time unit during the whole lifes- 
pan, the reproduction potential of age groups decreases also with age. In consequence, selection 
pressure on the oldest individuals will decrease which allows the accumulation of defects in the 
genes indispensable for surviving the later periods of the lifespan. Therefore, the accumulation 
of defective genes which are expressed after the minimum reproduction age should be observed 
and these defects would be responsible for the genetic program of ageing. 

In 1995 Penna described a simple computer model, in which he assumed that genes in 
genomes are switched on chronologically [2]. He observed that in the fraction of genes switched 
on before the reproduction age much less mutations were accumulated than in the fraction of 
genes switched on after the minimum reproduction age. It seems to be obvious that defects 
expressed before the reproduction kill the individual and they cannot be transferred to the 
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offspring. If defects are expressed after the minimum reproduction age, they can be inherited 
by the offspring and the chance of being transferred is still higher for genes expressed later 
during the lifespan. For description of the original model, its different versions and results of 
simulations of many demographic phenomena see [3], [1]. 

In the standard Penna model, the size of a population was controlled by Verhulst factor: 
V = 1 — N t /N max , where V describes the probability of survival for the individual independently 
of its age, N t corresponds to the actual size of the population and iV max is called the maximum 
capacity of the environment. It was the Verhulst factor which introduced the random death, 
assumed by Medawar, into the model. But in the other version of the model, the Verhulst 
factor determined the survival probability of the newborns only [3] in fact it controlled the 
birthrate. Like in nature, the reproduction probability and the survival of newborns depend 
on the resources availability and overcrowding. In simulations where the Verhulst test was 
set for newborns only, the random death of older individuals disappeared. In such models 
the accumulation of mutations in the fraction of genes expressed later during the lifespan is 
still observed. Furthermore, the genetic pool structure of such populations is better (lower 
fraction of defective genes). The results indicate that the random death during the lifespan in 
the population is not necessary for generating the nonrandom distribution of defective genes 
in the genomes, neither for emerging the genetically determined ageing. Using this model it 
is possible to recover in the computer simulations the age structure of populations which fits 
the real human populations [6]. Furthermore, it is possible to simulate the changes in the 
age structure of human populations and to predict demographic changes for the near future 
[7j, [S]. Nevertheless, such a modeling renders good results for fractions of populations in the 
reproduction ages, only. The age structure of the part of population before the reproduction 
age does not correspond to the real demographic data. 

The standard Penna model produces the low but constant fraction of defective genes in the 
set of genes expressed before the reproduction age. An important parameter of the model is 
the threshold number of defective phenotypes T. If this number is reached by an individual it 
has to die because of its genetic status. An individual cannot die before it reaches the threshold 
T e.g. for T = 3 the genetic death is impossible during the first two time units (Monte Carlo 
steps). Many authors tried to solve the problem of "the babies death" in the Penna model. 
They introduced lower threshold T for babies [5j, higher number of genes switched on in one 
time unit or introduced genes switched on in the periods between conception and birth [10J. 
All these modifications were biologically legitimated but there was no place in the model for 
more or less random relationships between individuals and environment. Such relations were 
introduced in the recent versions of the model pj], [12] where both, state of an individual and 
state of the environment fluctuate. If the sum of both fluctuations passes the limit set for 
homeodynamics of the organism then it has to die. If the babies are more sensitive for such a 
noise, their mortality is higher. 

This phenomenon wass very well seen in the ancient European human populations when 
many newborns and young children died because of infections. Recently, due to the intensive 
medical care, the mortality of newborns dropped. In the simulations described in this paper 
we have introduced some parameters mimicking the mechanisms which rescue the babies from 
the life threatening conditions. 
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2 The standard Penna model 



In the Penna diploid, sexual model each of N individuals of the population is represented by 
a genome composed of two bitstrings (haplotypes). Bits set for represent correct genes, bits 
set for 1 represent defective genes. Two bits at the same position in the bitstrings (locus) 
represent alleles. If both alleles are set for 1 then the phenotype determined by the locus is 
defective (i.e. defective genes are recessive). Loci are switched on chronologically. Organism 
dies if the declared number of defective loci (threshold T) has been switched on. The organism 
can reproduce if it survives until the age R (the reproduction age). During reproduction, each 
bitstring is copied and mutations are introduced into randomly chosen loci with a declared 
probability (0 is replaced by 1, bit set for 1 stays 1). Two new bitstrings recombine with a 
declared probability by exchanging their arms in the randomly chosen point. Each bitstring 
after these processes is a gamete. The offspring is produced by joining a gamete generated by a 
female with a gamete produced by a male. In the sexual model each newborn is declared to be 
female or male with the equal probability. The example of mortality curve and the distribution 
of defective genes generated in the standard Penna model simulations are shown in Fig. 1. 



3 The noisy Penna model 

For further extensions of the model we use only the diploid sexual Penna model described in 
the above section. In this standard model Verhulst factor controls the birthrate and there 
are no random deaths of organisms later during the lifespan. In the standard Penna model 
individuals die because of genetic death when they reach the threshold T of the expressed 
defective phenotypes, only. In the noisy version of the model there is no declared threshold 
T. Instead, we have introduced the fluctuations of the state of organisms. The energy of 
fluctuations increases with the number of switched on defective loci [llj. Thus the model is 

I i (t)=E(t) + P l (t), (1) 

where E(t) ~ N(nE{t), &1) corresponds to the fluctuations of environment in time t while 
Pi(t) ~ A^/xp^t), of (t)) corresponds to the inner fluctuations of individual i in time t. In the 
simplest case the expected value of both fluctuations is ^pjt) = ^E(t) = and the energy of 
fluctuations of the state of individual depends on its number of defective loci g^t) expressed 
till time t 

v 2 S)=°l+9M°l (2) 

Both models produce very similar results of simulations with characteristic gradient of 
defective genes expressed after the minimum reproduction age and very low mortality of the 
youngest individuals. The only difference between the two models concerns the mortality of the 
individuals during the first two time units. For threshold T = 3 in the standard model there 
are no genetic deaths during the first two time units (see Fig. 1). In the noisy model organisms 
may die even before the expression of any defect because of fluctuations (see Fig. 2). 
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Figure 1: The mortality curve (upper plot) and frequency of defective genes (lower plot) in 
the standard Penna model. Note the logarithmic y-axis of mortality plot, x-axis is scaled in 
the arbitrary age units. 
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Figure 2: The mortality curve (upper plot) and frequency of defective genes (lower plot) in the 
noisy Penna model. 



4 Mother care 

In the noisy Penna model the state of environment affects all individuals regardless of their 
age. We model the mother care as some kind of protection of the babies against the influence of 
fluctuations of the environment during the first periods of their lives. So, newborns are affected 
by states of the environment decreased by 

p(i,t) = 1 - exp(-(age(i,t) + 1)/A MC ). 

The individual dies if p(i,t)E(t) + P^t) > F. After \ MC steps the effect of mother care is 
negligible while in the early stages of life it is significant. 

We call this ,, mother care" to stress that this effect influences the very first periods of life 
but it could be the proper feeding of newborns with mother's milk as well as intensive medical 
care. In Fig. 3 we present results for \ M c = 4. The fraction of defects is similar for both 
models. 

5 Adaptation to the environmental conditions 

In the noisy Penna model the energy of fluctuations of individual is a sum of its inner noise 
and the noise of environment. The impact of these two components on the state and evolution 



5 



without mother care , 

with mother care 

10 20 30 40 50 60 

Figure 3: Mortality curves for the population with and without mother care. 



of population is different. The personal component is non-correlated and independent for each 
individual while the environmental noise is the same for each individual. That is the reason 
why the reactions of individuals for the environmental fluctuations are diversified. 

In the version of the noisy model presented above, the fluctuations have Gaussian distribu- 
tion with average /iE(t)®- Now, we introduce a signal into the expected state of environment. 
The signal /iE(t) is a periodical function with period D, thus /iE(t) — fi>E(t+D)- Individuals know 
that the signal is periodical, and are equipped with a mechanism of learning the signal. They 
estimate each component of the signal by weighted average of state of environment in survived 
periods. In the more formal way the learning mechanism affect the expected state of individual 
fluctuations 

oo 

^(0 = L ^ t - j * d)™^ -j* D ) 

3=1 

where L(i,t) = 1 if individual i lived at time t and otherwise while weights Wj are 

w j = e-W x -e- i ' x . 

This adaptation mechanism allows reducing the mortality in case when individual have 
learned the periodical signal. Results for different A are presented on Fig. 4. It is also ob- 
served in real populations, that mortality of newborns is higher than of a bit older individuals. 
The results depend on the maximal signal value (J,E(t) an d do not depend on the form of 
the periodic function, thus results for constant HE(t) — A are similar to those obtained with 
fJ>E(t) — ^4sin(t/7r) (results not shown). 

The next question is how the intensive protection of newborns against environment fluctua- 
tions could influence their mortality during the later periods of life. It is rather well known effect 
that children who are very strongly protected against any infections during the first periods of 
life and live in almost sterile conditions are more vulnerable for infections later. Plot shown in 
Fig. 6 indicate, that one could really expect slightly higher mortality of young individuals if 
they are isolated from environment influence during the very early periods of life. 
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Figure 6: The mortality curve for populations with and without mother care. For early periods 
of live the mortality of newborns is lower in case of mother care, but then it is a bit higher. 
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Figure 7: The mortality curve for the population with additional risk of death after H = 12 
years. 



6 Additional risk factors 

Suppose that after age H in every year an individual may die with some very small probability 
due to the random death (e.g. in a car accident or any other death of young people released 
from the parental care). In Fig. 7 we presented results for populations, in which individuals 
older than H = 12 dies in each year with probability p = 0.002 even if their state is lover than 
the F. The mortality curve resembles the curves observed in many real human populations. 
This random death introduced for individuals older than 12 could be also natural increase in 
mortality connected with reaching the puberty age. In Fig. 8 we presented the age structure of 
different human populations and for comparison, the mortality curves generated by simulations 
with properly rescaled the age axis. 

7 Conclusions 

The results of standard Penna model simulations reproduce the age distribution of the human 
population in its part after the minimum reproduction age. The noisy model enables simu- 
lation and analysis of parameters which influences the mortality of the youngest individuals. 
It will be possible to use this model for simulation the effect of pandemics and long distance 
immunological memory of populations. 
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Figure 8: The age structures (upper plots) and mortality curves of different human populations 
and populations generated in simulations with the extended noisy Penna model (lower plots). 



9 



References 



[1] P.B. Medawar, An unsolved problem of biology, An inaugural lecture delivered at University 
College London, published by H.K. Lewis&CO LTD London, (1952). 

[2] T.J. P. Penna, A bit-string model for biological aging. J. Stat. Phys. 78, 1629 (1995). 

[3] D. Stauffer, S. Moss de Oliveira, P.M.C. de Oliveira, J.S.Sa Martins, Biology, Sociology, Geology 
by Computational Physicists. Elsevier Amsterdam (2006). 

[4] A. Laszkiewicz, S. Cebrat, Monte Carlo Simulations of the Age Structure of the Human Popula- 
tion. J. Insurance Medicine 37(1), 3 (2005). 

[5] J.S. Sa Martins S. and Cebrat, Random deaths in a computational model for age-structured 
populations. Theory in Biosciences 119, 156, (2000). 

[6] E. Niewczas, S. Cebrat, D. Stauffer, The influence of the Medical Care on the Human Life 
Expectancy in 20th Century and the Penna Ageing Model. Theory in Biosciences, 119, 122-131, 
(2000). 

[7] A. Laszkiewicz, Sz. Szymczak, S. Cebrat, The Oldest Old and the Population Heterogeneity. Int 
J Mod Phys C, 14, 1355, (2003). 

[8] K. Bonkowska, S. Szymczak, S. Cebrat, Microscopic modeling the demographic changes. Int. J. 
Mod. Phys. C, 17(10), 1477, (2006). 

[9] K.N. Berntsen, Age-of-parent dependent mutation rate and weak children in the Penna Model 
in biological ageing. Int. J. Mod. Phys. C 7, 731 (1996). 

[10] E. Niewczas, A. Kurdziel, S. Cebrat, Housekeeping genes and death genes in the Penna ageing 
model. Int. J. Mod. Phys. C 11, 775 (2000). 

[11] P. Biecek, S. Cebrat, Fluctuations, Environment, Mutations Accumulation and Ageing. Int. J. 
Mod. Phys. C 17, 923 (2006). 

[12] P. Biecek, S. Cebrat. Immunity in the noisy Penna model. Int. J. Modern Phys. C. (in press). 



10 



